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Abstract 

We have carried out a new direct search for the CP violating decay Ks Sir^ 
with 1.7 fb~^ of e+e~ coUisions collected by the KLOE detector at the ^-factory 
DA$NE. We have searched for this decay in a sample of about 5.9 x 10^ KsKl events 
tagging the Ks by means of the Kl interaction in the calorimeter and requiring 
six prompt photons. With respect to our previous search, the analysis has been 
improved by increasing of a factor four the tagged sample and by a more effective 
background rejection of fake Kg tags and spurious clusters. We find no candidates 
in data and simulated background samples, while we expect 0.12 standard model 
events. Normalizing to the number of Ks — > 27r^ events in the same sample, we set 
the upper limit on BR{Ks -5- 37r°) < 2.6 x 10"^ at 90% C.L., five times lower than 
the previous limit. We also set the upper limit on the r/ooo parameter, |77ooo| < 0.0088 
at 90% C.L., improving by a factor two the latest direct measurement. 

Key words: e+e" collisions, DA$NE, KLOE, rare Ks decays, CP, CPT 



* Corresponding author 

Email addresses: stefano.iiiiscetti@lnf.infn.it (S. Miscetti), 
michal . silarskiOuj .edu.pl (M. Silarski). 



2 



1 Introduction 



The decay Kg — )■ 37r° violates CP invariance and its observation would be 
the first example of CP violation in Ks decays. The parameter r/ooo? the ratio 
of Ks to Kl decay amplitudes, is defined as: r^ooo = A{Ks — )• 37^^) /A{Kl — )■ 
37r°) = e+eQoo, where e indicates the Ks CP impurity and egoo the contribution 
of a direct CP- violating term. Since we expect egoo ^ ^ llJ; follows that 
r^ooo ~ e- In the Standard Model, therefore, BR{Ks 37r°) ~ 1.9 x 10"^ 
to a great degree of accuracy. The observation of such decay remains quite a 
challenge. 

Previous searches follow two alternative methods: via a fit to the interfer- 
ence pattern or via a direct search. The NA48 collaboration [2] has fit the 
Ks/Ki — i- 37r° interference pattern at small decay times finding 3? (t^ooo) = 
-0.002 ± O.Ollstat ± O.OlSsys and ^^(r^ooo) = -0.003 ± 0.013stat ± 0.017sys, 
corresponding to a limit on BR{Ks 37r°) < 7.4 x lO""^ at 90% C.L. The best 
upper limit on BR{Ks — > 3tc^) comes from the direct search performed by our 
experiment [3] based on 450 pb~^ of collision data collected during 2001-2002. 
We observed 2 candidates with an estimated background of 3.1 events in a 
pure tagged sample of 127 millions of Ks decays, that allowed us to quote a 
limit on BR{Ks 37r°) < 1.2 x 10"^ at 90% C.L. In this letter, we present a 
twofold improvement of this search based on a four times larger, and indepen- 
dent, data sample collected in 2004-2005 and on ameliorated techniques used 
for background rejection. 



2 The KLOE detector 



The KLOE experiment operated from 2000 to 2006 at DA<I>NE , the Frascati 
0-factory. DA<I>NE is an e^e~ collider running at a center-of-mass energy of 
~ 1020 MeV, the mass of the meson. Equal energy positron and electron 
beams collide at an angle of 7r-25 mrad, producing mesons nearly at rest. The 
detector consists of a large cylindrical Drift Chamber (DC), surrounded by a 
lead scintillating fiber Electromagnetic Calorimeter (EMC) both immersed in 
an axial 0.52 T magnetic field produced by a superconducting coil around the 
EMC. At the beams interaction point, IP, the spherical beam pipe of 10 cm 
radius is made of a Beryllium- Aluminum alloy of 0.5 mm thickness. Low beta 
quadrupoles are located inside the detector at a distance of about ± 50 cm 
from the interaction region. The drift chamber [1], 4 m in diameter and 3.3 m 
long, has 12582 all stereo drift cells with tungsten sense wires and is a really 
light structure, having the chamber shell made of carbon fiber-epoxy composite 
with an internal wall of ~ 1 mm thickness, and filled with a gas mixture of 90% 
helium, 10% isobutane, to minimize Ks regeneration and photon conversion. 
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The spatial resolutions are axy ~ 150 /im and cr^ ~ 2 mm. The momentum 
resolution is ct{p±)/p± ~ 0.4%. The calorimeter [5] covers 98% of the solid 
angle and is composed by a barrel and two endcaps, for a total of 88 modules. 
Each module is read out at both ends by photomultipliers for a total of 2440 
cells arranged in five layers. The energy deposits are obtained from the signal 
amplitude, while the arrival times and particles impact points are obtained 
from the spatial coordinates of the fired cell and the time differences. Cells close 
in time and space are grouped into energy clusters. The cluster energy E is 
calculated as the sum of the cell energies, while the cluster time T and position 
R are energy weighted averages. Energy and time resolutions are parametrized 
as aE/E = 5.7%/^E (GeV) and at = 57 ps/^E (GeV)©100 ps, respectively. 
The trigger |l6j uses both calorimeter and chamber information. In this analysis 
events are selected with the calorimeter trigger, requiring two energy deposits 
with £' > 50 MeV for the barrel and E > 150 MeV for the endcaps. Data are 
then analyzed by an event classification filter [7j, which selects and streams 
various categories of events in different output files. 

In this letter, we refer only to data collected during 2004-2005 for an integrated 
luminosity C = 1.7 fb~^ with the most stable running conditions and the 
best peak luminosity. A total of 5.1 billion mesons were produced, yielding 
1.7 X 10^ KsKl pairs. Assuming BR{Ks 37r°) ~ 1.9 x 10"^ about 3 signal 
events are expected to have been produced. 



3 Event selection 

At DA$NE the mean decay length of Kl, Al, is equal to ~ 340 cm and 
about 50% of Kls reach the calorimeter before decaying. A very clean Ks 
tag is provided by the interaction in the calorimeter (i^^-crash), which 
is identified by a cluster with polar angle 40° < 6cr < 140°, not associated 
to any track, with energy E^r > 100 MeV and with a time corresponding to 
a Kl velocity (3* ~ 0.2 in the rest frame. The average value of the e~^e~ 
center of mass energy W is obtained with a precision of 20 keV for each 200 
nb~^ running period using large angle Bhabha scattering events [8]. The value 
of W and the i^^-crash cluster position allows us to obtain, for each event, 
the direction of the Ks with an angular resolution of 1° and a momentum 
resolution of about 2 MeV. 

Because of its short decay length, Xs ~ 0.6 cm, the displacement of the Ks 
from the decay position is negligible. We therefore identify as photons from 
Ks decay, neutral particles that travel with (3 = 1 from the interaction point 
to the EMC ("prompt photons"). In order to retain a large control sample 
for the background while preserving high efficiency for the signal, we keep all 
photons satisfying E^ > 7 MeV and | cos 6*1 < 0.915. Each cluster is required 
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Eor(MeV) p 

Fig. 1. Distributions of the Kl energy deposit in the EMC {Ecr) and velocity in the 
(j) center of mass frame (/3*) for all events in the six photon sample. Black points 
represent data, while the MC backgromid simulation is shown as red histogram. 
The same distributions for events rejected by the track veto are shown in the inset 
and by the black triangles (data) and green filled histograms (MC simulation). 

to satisfy the condition — R^/c\ < min(3.5crf, 2 ns), where is the photon 
flight time and R the path length; at also includes a contribution from the 
finite bunch length (2-3 cm), which introduces a dispersion in the collision 
time. The photon detection efficiency of the calorimeter amounts to about 
90% for = 20 MeV, and reaches 100% above 70 MeV. After tagging the 
signal sample is selected requiring 6 prompt photons. For normalization we 
use the Ks — > 27[^ decay which is selected requiring 4 prompt photons. 
For both channels the expected background as well as the detector acceptance 
and the analysis efficiency are estimated using the Monte Carlo simulation 
tool of the experiment [7j. The simulation incorporates a detailed geometry 
and material composition of the KLOE apparatus and most of the data taking 
conditions of the experiment e.g. DA$NE background rates, position of the 
interaction point and beam parameters. All the processes contributing to the 
background were simulated with statistics twice larger than the data sample. 
Moreover, for the acceptance and the analysis efficiency evaluation a dedicated 
Ks — Bvr^ signal simulation was performed, based on a branching ratio equal 
to the best known upper limit [5] increased by a factor of 30 (about 5000 
events). 



3. 1 The six photon sample 

After Ks tagging, the selection of the Ks — ?• 37r'' decay is performed by 
searching six prompt photons from the decay of the pions. Requiring E^r > 
100 MeV, 0.17 < 13* < 0.28 and 40° < O^r < 140°, we count 76689 events. For 
these events we perform further discriminant analysis to increase the signal to 
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Fig. 2. Distribution of Xfu the tagged six photon sample for data ( black points), 
background simulation (solid histogram), and simulated Ks — )• Stt" signal (dashed 
histogram) . 

background ratio. 



The first analysis step aims to reject fake Ks tags. The distributions of Ecr 
and (3* for the selected data sample and background simulations are shown in 
Fig. [U In the /3* distribution, the peak around 0.215 corresponds to genuine 
Kl interaction in the calorimeter, while the flat distribution mainly originates 
from — )• KsKl — )• {Ks — )• Tr+Tr",/^^ — )■ 37r°) background events. In this 
case one of the low momentum charged pions spirals in the forward direction 
and interacts in the low-/3 insertion quadrupoles. This interaction produces 
neutral particles which simulate the signal of K^ interaction in the calorimeter 
(fake -ft^L-crash), while the K^ meson decays close enough to the interaction 
point to produce six prompt photons. To suppress fake K^-crash we first reject 
events having charged particles produced close to the interaction region (track 
veto). The distributions of the kinematical variables for the vetoed background 
events are shown in the inset of Fig. [H Taking advantage of the differences 
in the {3* and Ecr distributions between the tagged Ks events and the fake 
X^-crash, we have tightened the cuts on these variables: E^r > 150 MeV and 
0.20 < (3* < 0.225 {Ki-cTash hard). This improves by a factor 12 the rejection 
of this background with respect to the previous analysis [S]. 
The second source of background originates from wrongly reconstructed Ks — )■ 
27r° decays. The four photons from this decay can be reconstructed as six due 
to fragmentation of the electromagnetic showers (splitting). These events are 
characterized by one or two low-energy clusters reconstructed very close to 
the position of the genuine photon interaction in the calorimeter. Additional 
clusters can come from the accidental time coincidence between the decay 
and machine background photons from DA$NE. After tagging with the Kl- 
crash hard algorithm and applying the track veto we remain with a sample 
of about 50000 six photon events. A kinematic fit with 11 constraints has 
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Fig. 3. Distributions of events in the C37r-C27r plane, for six-photon sample tagged 
by KL-crash for data (left), and for the simulated Ks — >■ Sir^ decays (right). The 
boundaries of the background control regions Bl, B2, B3, B4, B5 and the signal 
region S are as specified in the text. 

been performed imposing energy and momentum conservation, the kaon mass 
and the velocity of the six photons in the final state. The distribution of 
the fit for data and background simulation, Xfu^ is shown in Fig. |2] together 
with the expected distribution for signal events. Cutting on xju considerably 
reduces the background from bad quality reconstructed events while keeping 
large signal efficiency. 

In order to improve rejection of events with split and accidental clusters, we 
have exploited the correlation between two x^-like variables named and 
Cstt- C2n is calculated by an algorithm selecting the best four out of six clusters 
satisfying the kinematic constraints of the two-body decay in the Ks — )■ — )• 
47 hypothesis. The pairing of clusters is based on the difference between the 
reconstructed 77 masses, rrii^^ and m2^^, and m^o P, and on the opening 
angle of the reconstructed pion directions in the Ks center of mass frame, 
^TT^, which is 180° for Ks — )■ 271^ events. Moreover, we check the consistency 
of the determination of the Ks four-momentum vector P^g determined from 
the reconstructed four-momentum of K^, with the sum of the photons four- 
momenta Prec = Y^i=i P'yt ■ For every possible pairing choice the algorithm 
calculates the (2-k defined as: 



_ (mi^^-m^o)^ (m2T,^-m^o) 

(,27r — n r 



a. 



2-K 
4 



a 



+ 



i=l 



CT, 



pIs-Ypi] [pks-i:/,.) 



i=l 



Px 



Vy 



Pz 



The minimization of (2tt gives the best two photon pairs fulfilling the Ks — )■ 
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SBOX 


Bl 


B2 


B3 


B4 


B5 


DATA 


220 ± 15 


5 ± 3 


15179 ± 123 


26491 ± 163 


6931 ± 83 


137 ± 12 


MC 


239 ± 11 


4 ± 3 


14905 ± 116 


26964 ± 169 


6797 ± 76 


100 ± 7 



Table 1 

Number of events populating control regions in the C37r-C27r plane defined in Fig. [3] 
after tight requirements on i^^^-crash and track veto. 

2tt^ — 7- 47 hypothesis. The resolutions used in Eq. [T] were estimated indepen- 
dently on data and MC simulation using a Ks — )■ 27r° — t- 47 control sample. 
The second x^-like variable, Csvr, instead verifies the signal hypothesis Ks — )■ 
371^ by looking at the reconstructed masses of the three pions. For each pair 
of clusters we evaluate Csvr as: 



a 



3tt 



a 



(2) 



As the best combination of cluster pairs, we take the configuration minimizing 
Cstt- The resolution on the 77 invariant mass in the 37r° hypothesis, a^^, was 
estimated applying the algorithm to the simulated Ks — )■ 37r° events. 
The distributions in the C37r-C27r plane for the data and Ks — )■ 37r° simulated 
signal are shown in Fig. [3l Signal events are characterized by small values of 
Cstt and relatively high (2tt- To compare data and Monte Carlo simulations we 
have subdivided the C37r-C27r plane into six regions Bl, B2, B3, B4, B5, and 
S as indicated in the left panel of FigJHl Region S, with the largest signal- 
to-background ratio, is the signal box, while B1-B5 are control regions used 
to check the reliability of the simulation and optimize our description of the 
experimental data. 

Simulation does not reproduce accurately the absolute number of events be- 
longing to different background categories. However, their kinematical prop- 
erties are reproduced quite well. To determine the background composition, 
and improve the description of experimental data, we have performed a binned 
likelihood fit of a linear combination of simulated C37r-C27r distributions to the 
same data distribution for all background categories [llj. The quality of the fit 
was controlled by comparing inclusive distributions of discriminating variables 
between data and simulation. Examples are presented in Fig. HI 
Table [T] shows the comparison of observed number of events with the expec- 
tations in each control region of the C37r-C27r plane. A satisfactory agreement is 
found. 

To further improve the Ks — )■ 2tt^ background rejection we cut on the dif- 
ference AE between the Ks energy determined from the reconstructed Ki 
four-momentum and the sum of energies E^- of the four prompt photons se- 
lected by the (2n algorithm: 



AE = K/2-5:E,J/as, 



(3) 
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Fig. 4. Inclusive distributions of the Cstt and C,2-k discriminating variables for six-pho- 
ton events: data (black points), background simulations (red curves). The dashed 
histograms represents simulated Ks 
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Fig. 5. Distributions of /S.E and Rmin discriminating variables for six-photon events: 
data (black points), background simulations (red curves). The dashed histograms 
represents simulated Ks — )• Svr*^ events. 

where o"^; is a 47 energy resolution estimated using the Ks — ?• 27r° — >■ 47 con- 
trol sample. For Ks — )■ 27r° decays accompanied by two background clusters, 
we expect AE ~ 0, while for Ks — )■ Stt*^ events AE ~ 111^^0/ a e- To further 
reject surviving Ks — J- 271" events with split clusters, we cut on the minimal 
distance between centroids of reconstructed clusters, Rmim considering that 
the distance between split clusters is on average smaller than the distance be- 
tween clusters originating from 7's of Ks — )■ 37r° decay. Distributions of these 
two discriminant variables are presented in Fig. [51 

Before opening the signal box, the cuts on the discriminant variables have 
been refined minimizing /c„t(x/it, (271, Cstt, AE, i?„,i„) = Nup/e^^, where €3^ 
stands for the signal selection efficiency and Nup is the mean upper limit 
(at 90% CL) on the expected number of signal events calculated on the 
basis of the corresponding expected number of background events B^xp = 
B exp{x}it yC2TT, Csn, AE , Rmin) from simulation ^llj. The outcome of the opti- 
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Til 
DO 




x3D 


Tl ATA 

_L/_rA. i Jr\. 


no _|_ in 


n -1- 1 


U ZIZ 1 


innn9-i- inn 




uy zc y 


MC 


93 ± 6 


± 0.06 


0.3 zb 0.3 


9862 ± 86 


6234 ± 71 


37 ± 4 


DATA 


13 ± 4 


± 1 


± 1 


1387 ± 37 


± 1 


± 1 


MC 


15 ± 2 


± 0.06 


± 0.06 


1194 ± 21 


± 0.06 


± 0.06 



Table 2 

Population of control boxes in the C37r-C27r plane after Xfu < 57.2 cut (first two 
rows) and before the last cut on Rmin (last two rows). 




50 100 150 200 250 300 

Fig. 6. Distributions of xju six-photon events in the signal box (left) and AE 
for six-photon events in the signal box applying the Xfu < 57.2 cut (right). Black 
points are data, background simulation is the red histogram. The dashed histogram 
represents simulated Ks — 37r'^ events. 



mizing procedure is xju < 57.2, AE > 1.88 and Rmin > 65 cm. The signal 
box is defined as: 4 < ^27r < 84.9 and Cstt < 5.2. Since the expected number of 
background events is estimated using the Monte Carlo simulations, we have 
checked at each stage of the analysis how well the simulation describes the 
experimental data. Distributions of Xfu: ^"^^ Rmin variables are presented 
in Fig. [6] and Fig. [7] for events in the signal box. In the right panel of Fig. [7| we 
present also the Rmin distribution just before the last cut Rmin > 65 cm. Ac- 
cording to the Monte Carlo simulation, these survived events are all Ks — )■ 
decays with two split clusters (95%), or one split and one accidental cluster 
(5%). Tab.[2]shows the number of observed (expected) events in the six regions 
of the {(2tt, Cstt) plane after the xju < 52.7 cut and before the cut on Rmin- A 
total efficiency of e^^ = 0.233 ± 0.012stat has been estimated. At the end of 
the analysis we find zero candidates in data and in the simulated background 
sample. 
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Fig. 7. Distributions of Rmin for six-photon events in the signal box applying the 
Xfit < 57.2 cut (left), and applying xju < 57.2 and AE > 1.88 cuts (right). Black 
points are data, background simulation is the red histogram. The dashed histogram 
represents simulated Ks Stt^ events. 

3.2 The normalization sample 



The Ks — )■ 27r normalization sample is selected requiring four prompt pho- 
tons. The Monte Carlo simulation shows a negligible amount of background of 
about 0.1% of the total. These events are essentially (p — K^K~ decays. After 
the /^L-crash hard tagging we find A^27r = (7.533 ± 0.018) x lO'' events. With 
the Monte Carlo simulations we have also determined the Kg — t- 27r'^ — )■ 47 
selection efficiency: €2-^ = 0.660 ± 0.002^^. The final number of produced 
Ks 27r° events is: Nnorm = N2^/e2n = (1-142 ± 0.005) x 10^. 



3.3 Evaluation of systematic uncertainties 



The systematic uncertainties for this search are related to the number of back- 
ground events and to the determination of the acceptance and selection effi- 
ciencies for the signal, esTr, and normalization, e2n, samples. 

For the tagged six-photon sample, we have investigated the uncertainties re- 
lated to the observed background at the end of the analysis. A difference of ~ 
2.4% in the EMC energy scale and resolution has been observed between data 
and MC simulation and has been studied using a control sample of Ks — )■ 271^ 
events. To evaluate the related systematic uncertainty on the background, we 
have repeated the upper limit evaluation with several values of the energy 
scale correction in the range of 2.2%-2.6%. Similarly, the analysis has been 
repeated modifying the resolution used in the definition of (271 and ^3,^ [TTj. 
Moreover, we have varied the AE and Rmin cuts by 5% and 6%, respectively. 
Cuts on Ecr and (3* were varied by ± 5% for both variables. The full analy- 
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Xfit cut 




1.46 


Rmin cut 




0.90 


TOTAL 


1.65 


2.30 



Table 3 

Summary table of the systematic uncertainties on the selection efficiencies for the 
signal, £37^, and normalization samples, e2-w 



sis was repeated in total twenty times applying each time one of the changes 
mentioned above. For all of these checks, we have observed no variation in the 
number of simulated background. 

For the acceptance of both the signal and normalization samples, we have 
evaluated the systematic uncertainty on the photon counting by comparing 
between data and simulation splitting, accidental probabilities and cluster re- 
construction efficiency. Another source of systematic uncertainties originates 
from the offline filter used, during data reconstruction, to reject cosmic rays 
and machine background events [llj. We consider instead negligible the influ- 
ence of trigger efflciency for both samples, since in our previous analysis |3] 
it was about 99.5% and the K^^-crash hard tagging requires a larger energy 
release in the calorimeter, which translates in an even larger trigger efflciency. 

The observed difference in the EMC energy scale and resolution between data 
and simulation enters also in the esjr evaluation. The effects have been esti- 
mated as Ae3,r/e37r = 1-0 % from the energy scale, and Ae3,r/e37r = 1.1 % 
from the resolution. The effect of the cut on Xj^^ has been tested construct- 
ing the ratio between the cumulative distributions for experimental data and 
simulation which leads to a systematics of Ae37r/e3^ = 1.46%. Finally, we have 
investigated the systematic effect related to the Rmin cut by varying its value 
by 6%, and estimated its contribution to be Ae3^/e37r = 0.9%. 

All the contributions to the systematic uncertainty are summarized in Tab. [3l 
with the total systematic uncertainty evaluated adding all effects in quadra- 
ture. 
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4 Results 



No events were observed on data in the signal region. Equally, no background 
events are found in the MC simulation based on twice the data statistics. In 
the conservative assumption of no background, we estimate an upper limit 
on the expected number of signal events \JL{Nev{Ks — )■ 37r°)) = 2.3 at 90% 
C.L., with a signal efficiency of e^^ = 0.233 ± 0.012^^4 ± O.OOGsys- In the same 
tagged sample we count Nnorm = (1-142 ± 0.005) x 10^ Ks — )■ 27r° events. 
Systematic uncertainties on background determination, as well as on the effi- 
ciency evaluation for the signal and normalization samples, are negligible in 
the calculation of the limit. 



Using the value BR{Ks 27r°) = 0.3069 ± 0.0005 [9] we obtain: 
BR{Ks 37r°) < 2.6 x 10"^ at 90% C.L. 



(4) 



which represents the best limit on this decay while improving by a factor of 
~ 5 our previous result 

This result can be translated into a limit on |?7oooh 



|'7ooo| 



A{Kl 37rO) 



tl BRjKs 37rO) 
\\~sBR{Kl ^37r^) 



< 0.0088 at 90% C.L. 



(5) 



This describes a circle of radius 0.0088 centered at zero in the 3f?(r7ooo)) ^(^7000) 
plane and represents a limit two times smaller than our previous result (Ref. [3]). 
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